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Abstract 

The now-famous Majumdar-Papapetrou exact solution of the Einstein- 
Maxwell equations describes, in general, A*' static, maximally charged 
black holes balanced under mutual gravitational and electrostatic inter- 
action. When N = 2, this solution defines the two-black-hole spacetime, 
and the relativistic two-center problem is the problem of geodesic motion 
on this static background. Contopoulos and a number of other workers 
have recently discovered through numerical experiments that in contrast 
with the Newtonian two-center problem, where the dynamics is completely 
integrable, relativistic null-geodesic motion on the two black-hole space- 
time exhibits chaotic behavior. Here I identify the geometric sources of 
this chaotic dynamics by first reducing the problem to that of geodesic 
motion on a negatively curved (Riemannian) surface. 
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1. The Majumdar-Papapetrou solution 



The general Reissner-Nordstrom metric 



for a charged black hole of mass M and charge Q takes a particularly simple 
form when the black hole is extremal with \Q\ = M: 



5 = - (1 - dt^ + (^1- dr^ + {dO^ + sinH dcp^) . (1) 

In the isotropic coordinates f = r — M, Eq. (1) can be written suggestively as 

l.f)"^,^+(l + ^)"™„ (2, 
where 

(3)jj = dr^+f^de^+sm^ed^^) 

denotes the flat Euclidean metric. The metric function 1 + M/f appearing in 
Eq. (2) has the form of a harmonic function in Euclidean space, and, mirac- 
ulously, when 1 + M/f is replaced with a more general harmonic function the 
metric Eq. (2) still remains a solution to the Einstein-Maxwell equations ([1-2]). 
More precisely, as first discovered by Majumdar and Papapetrou, the metric 

g = -U-^ dt^ + [dx'^ + dy^ + dz'^) (3) 

and the electromagnetic potential given by 

A = ±^dt (4) 

are a solution to the source-free Einstein-Maxwell equations as long as the func- 
tion U = U{x, y, z) satisfies Laplace's equation in flat space: 

3 

U,kk = ?7,XX + U,yy + U,^^ = . (5) 

fe=l 

Note that this solution is static {d/dt is a timelike Killing vector), but in general 
has no other symmetries. 

It was first realized by Hartle and Hawking ([3]) that with the choice 



i=i 1^ - 



for the potential U{x,y,z), the Majumdar-Papapetrou solution represents N 
extremal black holes, where the i'th black hole, stationary at the fixed position 

f = Ti, has mass Mi and charge \Qi\ = Mi. All charges Qi have the same 
sign given by the sign chosen in Eq. (4), which ensures that the holes remain in 
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equilibrium, balanced under mutual gravitational attraction and electrostatic 
repulsion. The apparent singularity in U{x, y, z) [and therefore in the metric 
Eq. (3)] at the positions f = fi is the usual coordinate singularity associated 
with static coordinates at an event horizon. Indeed, the surface area of a small 
coordinate sphere {t = const., |r — fi\ = e} around r = fi approaches the 
expected surface area of the horizon: 

hni (C/^47r|f - fi|^) = 47r hm. 

r—^fi r—^ri 

= 4:7rMi\ (7) 

and the metric can be extended analytically into the interiors of the black holes 
(into "negative |r — fi\") using Kruskal-like coordinates. As in the single black 

hole case (the extremal Reissner- Nordstrom solution) , the interiors of the black 
holes house true physical singularities where spacetime curvature blows up. 



(l + ^+0(l))^.-rl| 



2. Chaos in the two-black-hole spacetime 

When N = 2, the spacetime given by Eqs. (6) and (3) represents a rela- 
tivistic analogue to the two-center configuration in Newtonian gravity, in which 
the Newtonian gravitational field is generated by two point masses at fixed 
positions (i.e., the mutual gravitational interaction of the masses is ignored). 
Numerical investigations of null geodesic motion on this two-black-hole space- 
time by Contopoulos and coworkers ([4]) have revealed that the geodesies exhibit 
chaotic behavior in the vicinity of the two centers. More specifically, Contopou- 
los studies null geodesies whose spatial motion is confined to a two-dimensional 
symmetry plane; assuming the black holes are positioned along the z-axis, this 
plane is typically the surface {x = 0} (see Fig. 1). Numerical integration of the 
null geodesic equations then reveals that for geodesies that approach the black 
holes from infinity, it is essentially impossible to predict whether the orbit will 
plunge into the first hole, or the second one, or escape back out to infinity; in 
other words the qualitative behavior of the orbits near the black holes exhibits 
effectively stochastic features. This places the relativistic two-center motion in 
surprising contrast with the corresponding Newtonian problem (i.e., the motion 
of a massive test body in the gravitational field of two fixed centers) where the 
dynamics is known to be completely integrable (a classical result that goes back 
to the work of Jacobi and Liouville). 

In this paper I will argue that the chaotic behavior of the null geodesic 
flow has its roots in the spatial geometry of the two-black-hole spacetime, and 
I will do so by first showing that the dynamics of this flow can be reduced to 
that of ordinary geodesies on a negatively curved Riemannian surface. 
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3. Geometric analysis of the two-black-hole null geodesic flow 



I will rely on the well-known "Fermat's principle" in its relativistic for- 
mulation ([5]). Fermat's principle states that if M = RxT, is a static spacetime 
with metric 

g = goodf + ^^^h , 

where E is a 3-manifold, and goo < is a smooth function and ^^^h is a Ricman- 
nian metric on S (both independent of t) , then the null geodesies of ( Af . g) when 
projected onto S are precisely the Riemannian geodesies of the 3-geometry 

^^-^] ; (8) 



-goo 

and, furthermore, the afBne parameter (i.e., the arc length) along the projected 
geodesies in [S, ^^^h/{—goo)] is precisely the static time coordinate t measured 
along the null geodesies in {M, g) . In words that would have sounded familiar to 
Fermat, the principle states that light follows the path of shortest (or extremal) 
travel time between two given points in 3-space. 

In the multi-black- hole solution given by Eqs. (3) and (6), Fermat's prin- 
ciple shows that null geodesic flow in the asymptotically flat exterior region 
(outside the event horizons of the black holes) is equivalent to the Rieman- 
nian geodesic flow of the 3-geometry (S,/i), with the 3-manifold S given by 
Ti = \{N points}, and with the Riemannian metric /i on S given by 

h = = n'^idx^ + dy"^ + dz^) , (9) 

-goo 

where 

In the two-black-hole spacetime, I can assume without loss of generality that 
the holes are positioned along the z-axis at fi = (0,0,1) and r2 = (0)0,-1) 
(see Fig. 1). It is obvious that any two-plane containing the symmetry (z-) axis 
is a totally geodesic submanifold of E. As I will focus on null geodesies which 
lie (spatially) in such a symmetry plane, which I can assume to be the yz-plane 
{x = 0} as in Fig. 1, by Fermat's principle the null geodesic flow I need to study 
is equivalent to the geodesic flow on the two-dimensional Riemannian surface 
(S, h), where S = \ {(0, 1), (0, -1)}, and 

h = n^{dy^ + dz^), (11) 

n = fi+ , + , v. (12) 

Note that the geodesic flow of {S, h) corresponds, in the original spacetime, only 
to the null geodesic flow in the exterior of the black holes; null geodesic motion 
in the interior regions is not covered by this correspondence. This point will 
become clearer after a closer look at the topology and large-scale geometry of 
{S,h): 
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Global geometry of the Riemannian surface {S, h) 



Look closely at the behavior of the metric h near the centers, e.g., near 
r = fi = (0, 1). Introducing Euclidean polar coordinates {R, 9) centered around 
ri = (0, 1) (i.e., R = \ f— fi\), I can write the conformal factor O in the vicinity 
of fi as 



(13) 



and similarly I can write 



h = (^1 + ^ + ^ + 0{R)^ ' {dR' + R' dO^) . (14) 
Now introduce a new radial coordinate p = Mx jR. Then Eq. (14) becomes 

i + o(-) {dp^ + p^de^). (15) 
p ) 

A similar analysis can be carried out in the vicinity of the other center r = r-i 

with the same conclusion, namely that what looks like a singularity at r = ri 
(and similarly near the other center) is in fact an entire asymptotically flat 
Euclidean region squeezed into a small neighborhood of the "point" t\ = (0, 1) 
in the coordinate system {y,z). The global geometry of the surface (S, /i) is 
then as depicted in Fig. 2 below, with three asymptotically flat regions, one at 
r — > CX3, and two others at each of the centers f* — > ri and r ^ r-^. As a corollary, 
the surface (S', h) is geodesically complete. This is expected, since by Fermat's 
principle the afhne parameter (i.e., arc length) along the geodesies of {S,K) is 
the static time coordinate t measured along the null geodesies of the two-black- 
hole spacetime, and static time diverges to infinity at the event horizons of the 
black holes. In other words, a null geodesic in the two-black-hole spacetime falls 
into the i'th black hole if and only if the corresponding Riemannian geodesic in 
(S, K) escapes into the asymptotic region f ^ Vi. 

Local geometry of the Riemannian surface [S, h) 

The intrinsic geometry of a two-dimensional Riemannian manifold in the 
small is determined completely by the Gaussian curvature K (which is one-half 
the scalar curvature R). With the metric written in the conformally flat form 
Eq. (11), K is given by 

K = -i^A(logn) (16) 
= Jj[f23((Q-i)_^^ + (o-i)_^^)_n/_o_/] , (17) 

where A denotes the scalar Laplacian in the flat metric dy^ + dz^. It is straight- 
forward to compute K for the surface {S,h) by simply substituting Vt from 

Eq. (12) in Eq. (17). The result is a complicated expression, not particularly 
illuminating in its analytic form (which therefore I will not bother to give). A 
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plot of the curvature K as a. function of the coordinates y, z is given in Fig. 3 
(where I chose unit masses Mi = M2 = 1). It is apparent that K is strictly 
negative throughout S (and this is true for all masses Mi, M2 > 0). Both 
far away from and near the centers (where geometry is asymptotically flat) K 
approaches zero from below as expected (see Fig. 2). 



4. Can chaos in the two-black-hole spacetime be explained solely by 
the negatively curved geometry of {S, h)? 

In a Ricmannian manifold of arbitrary dimension n, negative sectional 
curvature causes neighboring geodesies to diverge exponentially ([6,7]). Recall 
the derivation of this well-known result: If Z denotes a vector field along the 
geodesic 7, Lie transported by a congruence of neighboring geodesies, then 

V^.V^.Z = i?^.2 7* , (18) 

where Rxy denotes the curvature operator V^Vy — VyVx — W[x,y]- Along 
7, an infinitesimal neighboring geodesic can then be defined abstractly as any 
solution of the "Jacobi equation;" a differential equation along 7 derived from 
Eq. (18). In a parallel-propagated basis {Eh} along 7 such that En =7*, Ja- 
cobi 's equation is 

d^Z"" 

~ —R'^nbnZ , (19) 

where s is the affine parameter. For a two-dimensional Riemannian manifold 
with Gaussian (= sectional) curvature K = J?i2i2, and with Z = Z^, Eq. (19) 
becomes 

£^ 

ds 

Assuming X < 0, and assuming the affine parameter s is small compared to 
\K / [dK / ds)\s=o\, Eq. (20) has generic solutions of the form 



2 - -KZ . (20) 




Z{s) ~ A(s)exp / ^/^ds') -FS(s)exp - / ^/^ds'] , (21) 



where A{s) and B{s) are slowly-varying amplitudes. It is clear that negative 
Gaussian curvature K results in an exponentially diverging Z{s) in general. 
When K{s) is bounded from above by a negative number, Eqs. (20)-(21) imply 
that 7 (as an orbit in the geodesic flow) has positive Liapunov exponents. Obvi- 
ously, exponential instability of orbits and positive Liapunov exponents are suffi- 
cient conditions for the presence of "sensitive dependence on initial conditions," 
the key ingredient of chaos. But are these criteria sufficient to demonstrate that 
chaotic behavior is indeed present? 

To investigate this question, let me briefly consider two examples from 
Newtonian gravitation. Recall that in classical mechanics, for a Hamiltonian 
system with Lagrangian function 

3,k 
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motion on a constant-energy (= constant-Hamiltonian) surface {H — E} is 
equivalent to geodesic motion on a Riemannian manifold, namely on the sub- 
manifold {V{q^) < E} of configuration space equipped with the Riemannian 
metric 



(Hamilton- Jacobi-Maupertius— • • principle; see [6]). Accordingly, in mathe- 
matical analogy with relativistic gravitation, so also in Newtonian gravity test- 
particle dynamics has a geometric description in terms of geodesic motion. In 
particular, motion in the Kepler and Newtonian two-center problems can both 
be described in terms of geodesies on a two-dimensional Riemannian surface, 
and this description can be put in exactly the same form as in Eqs. (11)-(12), 
except, of course, in the Kepler case the conformal factor Cl takes the form 



Plots of the Gaussian curvature of the metric [Eq. (22)] for the (planar) 
Kepler and Newtonian two-center problems are shown in Fig. 4. In both plots, 
E is chosen to be = —0.1 {E is chosen negative so that the geodesic flow 
describes the motion of bound orbits), and the masses are M = Mi = M2 = 1. 
There arc no surprises: As both systems are completely integrable with stable 
closed orbits, one would not expect negative curvature to be the dominant 
geometric feature. Indeed, in the Kepler case curvature is strictly positive, 
and in the two-center case it is mostly positive, with a small neighborhood of 
negative curvature in the vicinity of the centers; this small region of negative K 
corresponds to directional instabilities the orbits have while passing in between 
the two centers of attraction. [Note that in contrast with the black-hole surface, 
the center(s) in the Newtonian case are genuine singularities of the metric qe; 
however, these are not curvature singularities {K remains bounded as r — > fi), 
but rather conical singularities with a mass-independent angle deficit tt.] 

So far the association between negative Gaussian curvature on the one 
hand and chaotic behavior of the geodesic flow on the other appears to hold 
within the context of the three examples I discussed. Consider, however, one 
more example, this time the geodesic flow on the Riemannian surface S with 
only one extremal (Reissner-Nordstrom) black-hole; in other words with metric 
h given by Eq. (11) where O = (1 -|- Mjrf' ■ The Gaussian curvature K of the 
resulting geometry is plotted in Fig. 5 (with M = 1). As in the two-black-hole 
case (Fig. 3) , ii' is strictly negative everywhere. But the geodesic flow on this 
surface is a completely integrable system (angular momentum provides the sec- 
ond integral of motion). Clearly, then, negative curvature (sensitive dependence 
on initial conditions) is not sufficient for chaos: In fact, the unique closed (unsta- 
ble) geodesic in the geometry of Fig. 5 has strictly positive Liapunov exponents 




(22) 




(23) 



and in the Newtonian two-center problem it has the form 




(24) 
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as an orbit in the flow, so even the presence of positive Liapunov exponents does 
not always imply chaotic behavior. 

As others have done before, I would like to argue in this paper against the 
widespread practice in the physics literature of identifying chaos with merely 
the presence of positive Liapunov exponents. This is especially important in 
relativity (where there is no canonical choice for dynamical time) since whether 
or not a Liapunov exponent is positive depends crucially on the nature of the 
time parameter used in defining the exponent. In the next section I will describe 
a precise formulation for "chaos" [due to S. Willard ([8])] which I believe is 
particularly useful in relativity since it does not depend sensitively on the choice 
of time. In the following section (Sect. 6), I will demonstrate that null geodesic 
flow in the two-black-hole spacetime is chaotic according to this formulation. 



5. A precise formulation of chaos 

Central to our intuitive understanding of chaotic behavior is the notion of 

"sensitive dependence on initial conditions:" long-time prediction of motion in 
the phase space of a chaotic system is impossible since small initial perturbations 
of the orbits grow arbitrarily large as the system evolves in time. This, of course, 
is a vague idea in need of a precise mathematical formulation, and there exist 
various such formulations, the concept of Liapunov exponents being one of them. 
However, the exact content of our intuitive notion of sensitive dependence is not 
fully captured by the more precise concept of positive Liapunov exponents. For 
example, the phase flow {x = x, x G il"} has positive Liapunov exponents along 
all its orbits, but, clearly, this is not a chaotic system, and more complicated 
"counterexamples" with positive exponents can be found in which to discern 
that motion is non-chaotic would not be so easy. In order to conclude, on the 
basis of the presence of positive Liapunov exponents, that chaos is present, it 
is apparently necessary to make sure that the divergence of nearby orbits does 
not occur simply because these orbits escape to "infinity" under time evolution. 
What is needed to address this point is a mathematical formulation slightly 
more sophisticated than the concept of Liapunov exponents. 

Here, then, is my favorite "definition" of chaos, adopted from [8]: Restrict 
attention, for definiteness, to phase spaces M with metrizable topology. A 
dynamical system {M, (ft) is chaotic if it contains a "chaotic invariant subset," 
that is, a subset A c M such that: 

(CI) A is compact, and invariant under (ft, i.e., <Pt(A) C A yt G R. 
(C2) A has sensitive dependence on initial conditions. 
(C3) A is topologically transitive. 

The precise meaning of condition C2 (sensitive dependence on initial conditions) 
is the following: Fix a distance function p on M compatible with A^Ts topology. 
Condition C2 holds if there exists a fixed 6 > such that for all x € A and for 
every neighborhood U C A of x open in A, a point y & U and a, t > can be 
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found such that 

p[ipt{x),ipt{y)] > S . 

In other words, given any point a; G A, no matter how small a neighborhood 
U oi X I choose I can always find points y £ U O A whose orbits eventually 
diverge away from that of x under the flow ipt- Since A is compact, this notion 
of sensitive dependence on initial conditions is independent of the choice of p. 
Topological transitivity of A (condition C3) means the following: for every open 
U, V C A there exists ate R such that ^ptiU) fl V 5^ 0. 

Because the problem I study in this paper involves chaos "localized" in 
a bounded region of an asymptotically flat geometry (i.e., in the vicinity of 
the black holes), I will need to use a slightly generalized version of the above 
definition; my generalization is designed to be adapted to the essentially time- 
asymmetric nature of the problem (i.e., null geodesies approaching the black- 
hole region from infinity and plunging into the holes after exhibiting chaotic 
behavior). Namely, call a subset A C Ma, "chaotic future-invariant set" if 

(FCl) A is compact, and future-invariant under ^pt, i.e., i^t(A) C A Vt > 0. 

(FC2) A has sensitive dependence on initial conditions (defined as before). 

(FC3) A is topologically future-transitive. 

Note that topological transitivity of A as defined above (condition C3) ensures 
essentially that the flow is topologically "mixing;" this condition is designed 
to rule out situations in which A can be decomposed into multiple compact 
invariant sets. Clearly, topological transitivity would be an inappropriately 
strong condition to impose on a subset which is only future-invariant. Therefore, 
I modify this condition so as to demand that the flow on A is mixing only in the 
future direction, more precisely, I define A to be topologically future-transitive if 
there exists a time T > such that for every pair of open subsets U, V d ipriA) 
times t > Q and s > can be found such that ipt{U) fl ifis{V) ^- Clearly, a 
chaotic invariant set is also trivially a chaotic future-invariant set. The definition 
of a chaotic system can now be generalized to include any dynamical system 
which contains a chaotic future-invariant subset. 

Notice that this definition for chaos makes no reference to Liapunov ex- 
ponents; in fact, the rate of divergence of nearby orbits is not constrained in 
any way by the precise notion of sensitive dependence on initial conditions. 
This fact makes the definition especially interesting for applications in General 
Relativity: sensitive dependence as defined above holds for one choice of time 
function if and only if it holds for any other, as long as two choices of time are 
always related by a monotone-increasing diffeomorphism from the real axis R 
onto R. Of course, in general a mathematical definition is useful only if it is 
the subject of theorems, and there do exist theorems which demonstrate that 
many of the usual properties of chaotic systems can be derived from the above 
conditions C1-C3 (or FC1-FC3); I will not discuss these results here, but direct 
the reader to the literature, especially as listed in [8]. Instead I will turn now to 
the demonstration that the geodesic flow on the two-surface {S, h) (which, as I 
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discussed in Sect. 3, is equivalent to the null geodesic flow of the two-black-hole 
spacctimc) is chaotic according to the formulation of chaos I just described. It 
is important to note here that other studies (see [9]"[10] and references therein) 
have carried out this demonstration by searching for various more direct signa- 
tures of chaos in the two-black-hole geodesic flow; for instance, the existence of 
hyperbolic cycles and transverse homoclinic orbits in this flow is discussed in 
[9] , and the presence of positive Liapunov exponents is explored in [10] . 



6. "Proof of chaos in the two-black-hole null geodesic flow 

The geodesic flow on the Riemannian surface {S, h) can be described as 
a Hamiltonian dynamical system, with phase space M = the unit cotangent 
bundle of S, i.e., 

M = Ti*S = {{x,p) G T*S I Ibll = h'^'paPb = 1} , (25) 

and with the Hamiltonian fimction H{x,p) = \h"'^PaPb- I will denote the 
geodesic flow on T^S by the usual symbol </?f I showed in Sect. 3 that {S,h) 
has strictly negative Gaussian curvature, and recalled in Sect. 4 that negative 
curvature causes exponential divergence of the orbits in the geodesic flow. Now, 
if the surface S were compact, I could then simply define my invariant set A to 
be the entire phase space M — T^S: so chosen, A is compact when 5' is, and 
because of the negatively curved geometry of {S, h) , A has sensitive dependence 
on initial conditions, i.e., satisfies condition C2 as formulated in the previous 
section. It is not difiicult to show also that A is topologically transitive un- 
der the geodesic flow; therefore, if S were compact, all conditions C1-C3 for a 
chaotic invariant subset would be satisfied by this simple choice of A, i.e., the 
entire phase space would be a chaotic invariant set. Indeed, it is well known 
that geodesic flows of compact manifolds with negative sectional curvature are 
chaotic. (These flows in fact satisfy every criteria ever invented for chaos: they 
have positive Liapunov exponents, positive entropy, are mixing, are K-flows, 
• • • See [7] for an extensive but readable analysis of this classical problem.) The 
noncompactness of the two-black-hole Riemannian surface {S, h) is then the 
main difficulty I need to overcome in demonstrating the existence of a chaotic 
(future-) invariant subset in the (noncompact) phase space T^S. 

I will now construct a closed subset A cT*S which I claim is a chaotic 

future-invariant set for the geodesic flow. That A is compact and future- 
invariant will be evident from its construction, however, I will not be able to 
prove that A satisfies conditions FC2 and FC3. To prove these conditions, it 
would be sufiicient to combine the negatively curved geometry of (S, h) with the 
intricate topological structure that A appears to have; however, I cannot prove 
that A indeed has this intricate structure. As is usually the case with stiidies of 
chaotic behavior, the evidence for this structure is exclusively numerical. Some 
of this numerical evidence I will present here, and more of it can be found in 
the literature, e.g., in [4] and [10]. 
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First define subsets Fi, r2 and F of the phase space T^S as follows: Fj 
is the set of all points in T^S which fall into the i'th black hole as f — > oo, i.e., 

Ti = {me T^S I f[iptim)] — > fj as t ^ 00} , (26) 

and r is the set of all points which escape to the asymptotically flat region 
r = 00 as f — > 00, i.e., 

r = {meT^S\f[<ft{m)] — ^ 00 as 00}. (27) 

Since these subsets consist of points {x,p) such that the geodesic starting at x 
with initial tangent vector p eventually escapes to one of the three asymptotically 
flat regions of {S,h) (see Fig. 2), it is clear that both the Fj and F are open 
subsets in T^S. Also (and this will be important below), it is clear that F, Fi 
and F2 are mutually disjoint subsets, i.e., 

F n Fi = Fi n F2 = . 

Now define A as the closed subset 

A = complement (F U Fi U F2) = (F U Fi U F2)' ; (28) 

A is the set of all points which do not escape to any asymptotic region as 
f — > 00, i.e., the set of all future-imprisoned (e.g., periodic or quasi-periodic) 

orbits of the geodesic flow. This is obviously a future-invariant subset, but it 
is not necessarily compact (unless all imprisoned orbits are closed geodesies, 
which is not the case as numerical studies show). To cut A down to a compact 
size, introduce a compact subset D c S as follows (see Fig. 6): Draw a circle 
C in the asymptotically flat region r — > 00 which encloses both black holes, 
and draw circles Ci and C2 in the asymptotic regions r ^ fi and r ^ r2 
which enclose the black holes 1 and 2, respectively. Choose these circles large 
enough so that if n denotes the outward normal to C and Ci, a geodesic 7 which 
crosses any one of the circles in the outward direction [i.e., with /!,(7*,n) > 0] 
escapes to the corresponding infinity (and thus never crosses C or Ci again). 
By asymptotic flatness, it is clear that such circles C and Ci can be found (see 
Fig. 6). Now let D c S be the compact region bounded by the circles, in other 
words, define D to be the unique connected component of S \ (C U Ci U C2) 
such that = CU Ci U C2. Then put 

A = A n T^D . (29) 

So constructed, A is clearly both compact (a closed subset of a compact set) 
and future-invariant. I claim that this A C S is a chaotic future-invariant 
subset for the geodesic flow on T*S. 

As I mentioned above, I am not able to prove that A satisfies the con- 
ditions FC2 and FC3 of Sect. 5. Nevertheless, a great deal of insight into the 

structure of A can be obtained by numerically integrating the geodesic equations 
on {S,h). Extensive numerical studies of this kind have been reported in [4] 
and [10]. Although I will base the following observations on my own minimal 
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investigation of the (numerical) structure of A, these observations are supported 
by the more extensive numerical evidence already published in the literature. 

Because of the exponential instability of all orbits in the geodesic flow, 
it is clear that a direct computer proof of the existence of a future-imprisoned 
orbit (lying in A) is impossible: any real orbit in the computer will eventually 
diverge away from A because of numerical instabilities, even if initially it lies 
in A. So it might appear at first that by relying on numerical integration it is 
impossible to even prove that A is nonempty! This is not the case, however; 
numerical integration does yield an indirect proof that orbits which lie in A 
exist. More precisely, consider those orbits whose starting points are on the 
z-axis and whose initial (unit) tangent vectors are entirely in the ^/-direction. 
[See Fig. 7; all orbits plotted in Fig. 7 are of this kind. Also, although the orbits 
plotted in Fig. 7 are (mostly) with unit masses Mi = M2 = 1 and (some) with 
masses Mi = 2, M2 = 1, similar behavior is observed with all positive choices 
of Ml, M2.] In the following, I will not make any distinction between points 
on the 2;-axis and initial conditions for the orbits in S; the initial-tangent- 
vector part of the initial conditions is fixed throughout to be a unit vector in 
the y direction. Now, by numerically integrating these orbits into the future, 
the following features can be observed: (i) Consider any open interval of initial 
conditions (starting points) on the 2:-axis lying in the vicinity of the centers. 
No matter how small this interval is, there are always points in it which belong 
to r, Fi and (ii) In any such interval, between any two points that belong 
to a distinct pair of the subsets F, Fi and F2, there exists a third point which 
belongs to the subset other than the two in the pair. 

Note that since F, Fi and F2 are open sets, both of the statements (i) 
and (ii) are "stable" numerically, i.e., they can be verified with arbitrarily- high- 
accuracy numerical calculations. Already, the statement (i), combined with the 
observation that F and Fj are disjoint, proves that A is nonempty: a connected 
open interval in R cannot be the union of three disjoint open subsets, therefore, 
in any interval of the kind described in (i) there must exist points which belong 
to A. As I remarked above, to prove that A satisfies the conditions FC2 and FC3 
of Sect. 5, it is suflacient to combine the exponential instability of the geodesic 
flow on TiS with the everywhere-dense topological structure of A, i.e., the 
structure of a Cantor set of periodic or quasi-periodic orbits, so that every open 
neighborhood of any point m G A contains points of A other than m. That 
A indeed has this structure is strongly suggested by the numerical evidence 
discussed here and more extensively in [4] and [10]. However, the discovery of 
an analytical proof of this topological structure remains an open problem. 
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FIGURE CAPTIONS 



Figure 1. Contopoulos's (and also this paper's) analysis of the two-black- 
hole null-geodesic flow is confined to those null geodesies which lie in a two- 
dimensional surface of symmetry such as the yz-plane {x — 0}. 

Figure 2. The geometry of the Ricmannian manifold (S, h) in the large. Note 
that this is not the actual geometry of the surface {x — 0} in the physical metric 
g on the two-black-hole spacetime, but, rather, it is the physical geometry with 
an extra conformal factor introduced in accordance with Fermat's principle. 
In particular, only the asymptotic region r ^ oo corresponds to the usual 
asymptotic region in the physical spacetime; the asymptotic regions r — s- 
exist because of the singular behavior of the static time coordinate t at the 
event horizons of the black holes. Accordingly, a null geodesic in the two-black- 
hole spacetime falls into the I'th black hole if and only if the corresponding 
Riemannian geodesic in {S,h) escapes into the asymptotic region r — + r,. 

Figure 3. The Gaussian curvature K of the Riemannian surface [S, h) as a 

function of the coordinates (y, z). The masses are chosen to be Mi = M2 = 1 for 
this plot; but the qualitative features of K are identical for all positive masses. 
In particular, K is strictly negative throughout S, and approaches zero in all 
three asymptotic regions, i.e., both as r ^ 00 and as r — > ri, i = 1, 2. 

Figure 4. No surprises in Newtonian gravity for the connection between nega- 
tive curvature and chaotic geodesic motion: With completely integrable geodesic 
flows, Gaussian curvature of the metric qe [Eq. (22)] is positive for both the Ke- 
pler problem (plot on the left; strictly positive K) and the Newtonian two-center 
problem (plot on the right; K positive except in a small neighborhood of the 
centers). 

Figure 5. The Riemannian surface S with the metric corresponding to the null 
geodesic flow of a single extremal Reissner-Nordstrom black hole has strictly 
negative Gaussian curvature, and the Liapunov exponent of its (unique) closed 
geodesic is positive. But there is no trace of chaos here: with angular momentum 
as the second integral of motion, the geodesic flow of this surface is a completely 
integrable Hamiltonian system. 

Figure 6. Construction of the compact set D C S used in defining the compact 
future-invariant subset A C TiS [see Eq. (29)]. The circles C and Ci are chosen 
large enough so that any geodesic crossing them in the outward direction never 
comes back (it escapes to the corresponding asymptotic infinity). The subset D 
is the compact connected region bounded by the three circles. 

Figure 7. Closed (or almost closed) orbits in the geodesic flow on {S,h). The 
top four plots are drawn with unit masses Mi = M2 = 1, and the two plots 
at the bottom of the flgure are drawn with masses Mi = 2 and M2 = 1. As 

the orbits get more complicated (and therefore their periods become longer), 
numerical instabilities set in as soon as or before the full shape of the orbit 
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becomes apparent (as happens in the middle two plots). Recall that all these 
orbits are unstable because of the negatively curved geometry of {S, h). 
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